Time discretization and quantization methods for 
optimal multiple switching problem* 

Paul Gassiat 1 ), Idris Kharroubi 2 \ Huyen Pham 1 -'' 3 '' 

September 23, 2011 
Revised version: February 7, 2012 



^ Laboratoire de Probabilites et 

Modeles Aleatoires, CNRS, UMR 7599 

Universite Paris 7 Diderot, 

gassiat, pham at math.univ-paris-diderot.fr 
3 ) CREST-ENSAE 

and Institut Univcrsitairc de France 



2 ) CEREMADE, CNRS, UMR 7534 
Universite Paris Dauphinc 
kharroubi at ceremade.dauphhic.fr 



Abstract 

In this paper, we study probabilistic numerical methods based on optimal quanti- 
zation algorithms for computing the solution to optimal multiple switching problems 
with regime- dependent state process. We first consider a discrete-time approximation 
of the optimal switching problem, and analyze its rate of convergence. Given a time 
step h, the error is in general of order (/ilog(l//i)) 1 ^ 2 , and of order h 1 ^ 2 when the 
switching costs do not depend on the state process. We next propose quantization 
numerical schemes for the space discretization of the discrete-time Euler state process. 
A Markovian quantization approach relying on the optimal quantization of the normal 
distribution arising in the Euler scheme is analyzed. In the particular case of uncon- 
trolled state process, we describe an alternative marginal quantization method, which 
extends the recursive algorithm for optimal stopping problems as in [2J. A priori L p - 
error estimates are stated in terms of quantization errors. Finally, some numerical tests 
are performed for an optimal switching problem with two regimes. 
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1 Introduction 



On some filtered probability space (f^J 7 , F = (.Ft) t>0 ,P), let us introduce the controlled 
regime-switching diffusion in M. d governed by 

dX t = b{X t ,a t )dt + a{X t ,a t )dW t , 

where W is a standard d-dimensional Brownian motion, a = (r n , i n ) n E A is the switching 
control represented by a nondecreasing sequence of stopping times (r n ) together with a 
sequence (i n ) of F Tn -measurable random variables valued in a finite set {1, . . . , q}, and at 
is the current regime process, i.e. at = i n for r n <t < r n+ i. We then consider the optimal 
switching problem over a finite horizon: 



j-T 

Vo = sup E 



rT 

/ f(X t ,at)dt + g(X T ,a T )- V] c(X TfI , . (1.1 

■/ o _ ^ J 



r„<T 



Optimal switching problems can be seen as sequential optimal stopping problems belonging 
to the class of impulse control problems, and arise in many applied fields, for example in real 
option pricing in economics and finance. It has attracted a lot of interest during the past 
decades, and we refer to Chapter 5 in the book [T7] and the references therein for a survey 
of some applications and results in this topic. It is well-known that optimal switching 
problems are related via the dynamic programming approach to a system of variational 
inequalities with inter-connected obstacles in the form: 

min --r£- K x : i)-D x Vi - -tr(a(x, i)a(x, i)' D 2 x Vi) - f(x, i) , (1.2) 

Vi — max(vj — c(x,i,j)) = on [0, T)xl d , 

together with the terminal condition Vi(T,x) = g(x,i), for any i = l,...,q. Here Vi(t,x) 
is the value function to the optimal switching problem starting at time t € [0, T] from the 
state Xt = x E M rf and the regime at = i E {1, . . . , q}, and the solution to the system (jl.2p 
has to be understood in the weak sense, e.g. viscosity sense. 

The purpose of this paper is to solve numerically the optimal switching problem Ql.lj) . 
and consequently the system of variational inequalities (jl.2p . These equations can be solved 
by analytical methods (finite differences, finite elements, etc ...), see e.g. [14], but are known 
to require heavy computations, especially in high dimension. Alternatively, when the state 
process is uncontrolled, i.e. regime-independent, optimal switching problems are connected 
to multi-dimensional reflected Backward Stochastic Differential Equations (BSDEs) with 
oblique reflections, as shown in [9] and pH], and the recent paper [5] introduced a discretely 
obliquely reflected numerical scheme to solve such BSDEs. From a computational view- 
point, there are rather few papers dealing with numerical experiments for optimal switching 
problems. The special case of two regimes for switching problems can be reduced to the re- 
solution of a single BSDE with two reflecting barriers when considering the difference value 
process, and is exploited numerically in [8]. We mention also the paper [3], which solves an 
optimal switching problem with three regimes by considering a cascade of reflected BSDEs 
with one reflecting barrier derived from an iteration on the number of switches. 
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We propose probabilistic numerical methods based on dynamic programming and opti- 
mal quantization methods combined with a suitable time discretization procedure for com- 
puting the solution to optimal multiple switching problem. Quantization methods were 
introduced in [2] for solving variational inequality with given obstacle associated to optimal 
stopping problem of some diffusion process (Xt). The basic idea is the following. One first 
approximates the (continuous-time) optimal stopping problem by the Snell envelope for the 
Markov chain (Xt k ) defined as the Euler scheme of the (uncontrolled) diffusion X, and then 
spatially discretize each random vector Xt k by a random vector taking finite values through 
a quantization procedure. More precisely, (X tk )k is approximated by (Xk)k where Xk is 
the projection of Xt k on a finite grid in the state space following the closest neighbor rule. 
The induced /^-quantization error, \\X tk — Xk\\ p , depends only on the distribution of X tk 
and the grid, which may be chosen in order to minimize the quantization error. Such an 
optimal choice, called optimal quantization, is achieved by the competitive learning vector 
quantization algorithm (or Kohonen algorithm) developed in full details in [2]. One finally 
computes the approximation of the optimal stopping problem by a quantization tree algo- 
rithm, which mimics the backward dynamic programming of the Snell envelope. In this 
paper, we develop quantization methods to our general framework of optimal switching 
problem. With respect to standard optimal stopping problems, some new features arise 
on one hand from the regime-dependent state process, and on the other hand from the 
multiple switching times, and the discrete sum for the cumulated switching costs. 

We first study a time discretization of the optimal switching problem by considering 
an Euler-type scheme with step h = T/m for the regime-dependent state process (Xt) 
controlled by the switching strategy a: 

X tk+1 = X tk +b(Xt k ,at k )h + a(X tk ,a tk )Vh$k+i, t k = kh, k = 0, . . . ,m, (1.3) 

where k = 1, . . . ,m, are iid, and A/"(0, /^-distributed. "We then introduce the optimal 
switching problem for the discrete-time process (Xt k ) controlled by switching strategies 
with stopping times valued in the discrete time grid {t^,k = 0, . . . ,m}. The convergence 
of this discrete-time problem is analyzed, and we prove that the error is in general of order 
(hlog(l/h)) 2 , and of order ,as for optimal stopping problems, when the switching costs 
c(x,i,j) = c(i,j) do not depend on the state process. Arguments of the proof rely on 
a regularity result of the controlled diffusion with respect to the switching strategy, and 
moment estimates on the number of switches. This improves and extends the convergence 
rate result in [5] derived in the case where X is regime- independent. 

Next, we propose approximation schemes by quantization for computing explicitly the 
solution to the discrete-time optimal switching problem. Since the controlled Markov chain 
(X tk )k cannot be directly quantized as in standard optimal stopping problems, we adopt a 
Markovian quantization approach in the spirit of [15], by considering an optimal quantiza- 
tion of the Gaussian random vector fik+i arising in the Euler scheme (jl.3p . A quantization 
tree algorithm is then designed for computing the approximating value function, and we 
provide error estimates in terms of the quantization errors — ~&k\\p and state space grid 
parameters. Alternatively, in the case of regime-independent state process, we propose a 
quantization algorithm in the vein of [5] based on marginal quantization of the uncontrolled 
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Markov chain (Xt k )k- A priori X^-error estimates are also established in terms of quantiza- 
tion errors \\Xt k — Xk\\ p . Finally, some numerical tests on the two quantization algorithms 
are performed for an optimal switching problem with two regimes. 

The plan of this paper is organized as follows. Section 2 formulates the optimal swit- 
ching problem and sets the standing assumptions. We also show some preliminary results 
about moment estimates on the number of switches. We describe in Section 3 the time dis- 
cretization procedure, and study the rate of convergence of the discrete-time approximation 
for the optimal switching problem. Section 4 is devoted to the approximation schemes by 
quantization for the explicit computation of the value function to the discrete-time optimal 
switching problem, and to the error analysis. Finally, we illustrate our results with some 
numerical tests in Section 5. 

2 Optimal switching problem 
2.1 Formulation and assumptions 

We formulate the finite horizon multiple switching problem. Let us fix a finite time T 
£ (0,oo), and some filtered probability space (f^J 7 , F = (J c t) t>0> P) satisfying the usual 
conditions. Let I q = {1, ...,(/} be the set of all possible regimes (or activity modes). 
A switching control is a double sequence a = (r n ,i n ) n >o, where (r n ) is a nondecreasing 
sequence of stopping times, and i n are Jv n -measurable random variables valued in I q . The 
switching control a = (r n , i n ) is said to be admissible, and denoted by a £ A, if there exists 
an integer-valued random variable N with tn > T a.s. Given a = (r„, i n ) n >o £ A, we may 
then associate the indicator of the regime value defined at any time t £ [0, T] by 

It = i-ol{0<t<T } + X/ nl { r "<*< T ™+i}' 
n>0 

which we shall sometimes identify with the switching control a, and we introduce N(a) the 
(random) number of switches before T: 

N(a) = #{n > 1 : r n < T} . 

For a € A, we consider the controlled regime-switching diffusion process valued in M. d , 
governed by the dynamics 

dX s = b{X s ,I s )ds + a{X s ,I s )dW s , X = x eM d , (2.1) 

where W is a standard d-dimensional Brownian motion on (f^J 7 , F = (J~t)o<t<T, F)- We 
shall assume that the coefficients b { = b(.,i): R d R d , and cij(.) = a(.,i) : M d ->• R dxd , i 
£ I q , satisfy the usual Lipschitz conditions. 

We are given a running reward, terminal gain functions /, g : M d X I q — >■ M, and a cost 
function c : R d x I q x I q ->■ R, and we set = f(.,i), gi{.) = g(.,i), (Hj(.) = c(.,i,j), i,j 
£ I q . We shall assume the Lipschitz condition: 

(HI) The coefficients fi, gi and Cij, i,j £ I g are Lipschitz continuous on R d . 
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We also make the natural triangular condition on the functions Cj,- representing the 
instantaneous cost for switching from regime i to j: 

(He) 

ca(.) = 0, i G Ig, 
inf dj(x) > 0, for i,j G I q , j / i, 

x&R d 

inf [cij{x) + Cj k (x) - Ci k (x)] > 0, for k G I q , j ^ i, k. 



The triangular condition on the switching costs in (He) means that when one changes 
from regime i to some regime j, then it is not optimal to switch again immediately to 
another regime, since it would induce a higher total cost, and so one should stay for a while 
in the regime j. 

The expected total profit over [0, T] for running the system with the admissible switching 
control a = (r n , i n ) S A is given by: 



'0 n=l 
The maximal profit is then defined by 



J (o) = E[ / f(X t ,I t )dt + g(X T ,I T )- Vc(x Tn , tn _ 1)tn ; 



Vb = sup Jo (a). (2.2) 

The dynamic version of this optimal switching problem is formulated as follows. For (t, i) 
€ [0, T] x 1^, we denote by At,i the set of admissible switching controls a = (r n , i n ) starting 
from i at time t, i.e. to = Lq = i. Given a E At,i, and x G M. d , and under the Lipschitz 
conditions on b, a, there exists a unique strong solution to (|2.ip starting from x at time i, 
and denoted by {Xt' x ' a ,t < s <T}. It is then given by 

fT„+iAs [T n+ iAs 

X*™ = x+J2 b Ln (X t i mdu+ a tn (X^ a )dW u ,t<s<T. (2.3) 

The value function of the optimal switching problem is defined by 

T N{a) 

Vi(t,x) = supeT/ f(X^ a J s )ds + g(X^ a J T )-Y / <K^ a ^n-utn)}, (2.4) 
aeAt,i l Jt n=1 J 

for any (t,x,i) G [0, T] x W 1 x I g , so that Vb = maxj G i 9 Vi(0,xo). 
For simplicity, we shall also make the assumption 

9i(x) > max[gj(x) - Cij(x)], V(x, i) G M d x I„. (2.5) 

This means that any switching decision at horizon T induces a terminal profit, which is 
smaller than a no-decision at this time, and is thus suboptimal. Therefore, the terminal 
condition for the value function is given by: 



Vi(T,x) = gi(x), (x,i)€R d xI q 
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Otherwise, it is given in general by Vi(T,x) = maxj & i q [gj(x) — Cij(x)}. 

Notations. |.| will denote the canonical Euclidian norm on M. d , and (.|.) the corresponding 
inner product. For any p > 1, and Y random variable on (fi, J 7 , P), we denote by \\Y\\ p = 
(K\Y\ p )p. 

2.2 Preliminaries 

We first show that one can restrict the optimal switching problem to controls a with 
bounded moments of N(a). More precisely, let us associate to a strategy a G At,i, the 
cumulated cost process C t,x,a defined by 

C t,x, a = ^ c ( X £f>Vn-l,^)lr„<u, t<U<T. 
n>l 

We then consider for i£l d and K > the subset Af^x) of At,% defined by 

A£i(x) = {a E A,i :E|C^ x ' a | 2 < A'(l + |x| 2 )}. 

Proposition 2.1 Assume that (HI) and (He) holds. Then, there exists some positive 
constant K s.t. 



T N(a) 

Vi (t,x) = sup E / /(I^M^dHg^./r) - V c(I^V n _ lltn );2.6) 
for any (t,x,i) G [0,T] x R d x I q . 



Remark 2.1 Under the uniformly strict positive condition on the switching costs in (He), 
there exists some positive constant ?7 > s.t. N(a) < nC^' for any (t, x, i) G [0,T]xM d xI g , 
a G At,i- Thus, for any a G Af^x), we have 



E\N(a)\ 2 < r]K(l + 



x 



which means that in the value functions Vi(t,x) of optimal switching problems, one can 
restrict to controls a for which the second moment of N(a) is bounded by a constant 
depending on x. 

Before proving Proposition 12.11 we need the following Lemmata. 
Lemma 2.1 For all p > 1, there exists a positive constant K p such that 

< K p (l + \x\), 



sup 

for all (t, x, i) G [0, T] x R d x I 



sup \Xl' x ' a \ 
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Proof. Fix p > 1. Then, we have from the definition of xl' x ' a in (|2.3p . for (t,x,i) £ 
[0,T] xK d xI g ,a£ 



E 



sup \X l s ' x ' 
s£[t,r] 



t.x.a \P 



r fTn+lhr 



+ E 



sup Yl 

seM Tn<! 



r„<r T ™ 

r„ +1 As 



for all r £ [i, T]. From the linear growth conditions on and <7j, for i £l q , and Burkholder- 
Davis-Gundy's (BDG) inequality, we then get by Holder inequality when p > 2: 



E 



sup \xl' x ' a \ p 

L se[t,r] 



< K p [l + \x\ p + I E 



sup \xl' x ' a \ p du 

s£[t,u] 



for all r £ [t,T]. By applying Gronwall's Lemma, we obtain the required estimate for p > 
2 , and then also for p > 1 by Holder inequality. □ 



Lemma 2.2 Under (HI) and (He), the functions Vi, i £ I q , satisfy a linear growth con- 
dition, i.e. there exists a constant K such that 



i(t,x)\ < K(l + \x\), 



for all (t, x, i) £ [0, T] x R d x I 



Proof. Under the linear growth condition on fa, g^ in (HI), and the nonnegativity of the 
switching costs in (He), there exists some positive constant K s.t. 



E 



N(a) 



'-Jt 

< k(i + e 



fiXt^J^ds + giX^Jr) - £ cpf^'Vn-i,^ 



n=l 



t,x,a I 



sup | X^ 

u£[0,T] 

for all (t, x, i) £ [0, T] x R d x Iq, a £ ^4t,i. By combining with the estimate in Lemma |2,H 
this shows that 

Vi{t,x) < K(l + \x\). 
Moreover, by considering the strategy a with no intervention i.e. N(a°) = 0, we have 



Vi(t,x) > E I f(Xl> x > a ,i)ds + g(X l /' a ,i) 



> -K 1 + E 



t,x,a I 



sup \X% 

ue[o,T] 



Again, by the estimate in Lemma |2. 11 this proves that 

Vi(t,x) > -K(l + \x\) , 
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and therefore the required linear growth condition on V{. □ 
We now turn to the proof of the Proposition. 

Proof of Proposition 12.11 The proof is done in 4 steps. Given a € At.i, we will denote 

T N(a) 
J(t,x,i;a) = E / f(X t s ^°',I s )ds+g(X t /' a ,I T )-J2< X ^' a ^n-Uh 

n=l 

• Step 1. First, we notice that the supremum in the definition of Vi(t, x) may be taken over 
A s ti , where 

= I a = (T n ,i n ) £ A t ,i : (r n ) is strictly increasing }. 

Indeed, it is always suboptimal to switch several times at a single date due to the triangular 
condition (He). 

• Step 2. We now prove that it is enough to take the supremum over the strategies in A s t ^ , 
where 



AT = {*eAt ti :E\C^\ 2 < +oo}. 



For any a = (Tk,i-k)k>o G A^, define a n = (^\^)fc>o as the strategy obtained from a by 
only keeping the first n switches, i.e. 

tJ? = oo, k > n 

Note that for each n, a n £ A% . Now since a and a n (and the associated processes) 
coincide on {N(a) < n}, and by positivity of the switching costs, 



J(t, x, i; a) — J(t, x, i; a n ) 
< e' ' 



J (f(X^ a ,I s ) - f(Xt>*> an , I s ))ds + g(X^ a ,I T ) - g(X^ an , I T )) 1 
< K(l + \x\)F(N(a) > n) 1/2 , 



{N(a)>n} 



by Cauchy-Schwarz inequality, linear growth of /, g and Lemma l2.1i Hence letting n — > oo, 
and since N(a) < oo a.s., we obtain 

J(t, x,i;a) < lim inf Jit, x, i; a n ), 

n— too 

which proves the required assertion. 

• Step 3. To each a E A 6 t '°°, we associate the process (Y t ' x,a , Z t,x,a ) solution to the following 
Backward Stochastic Differential Equation (BSDE) 

= g(X t /> a ,I%)+ F f(X^ a ,I«)ds (2.7) 

J u 

T 

Z t,x,a dWs _ C hx,a + c t,x,a _ t <U<T 
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and satisfying the condition 



E 



sup | Y t 

se[t,T] 



t,x,a\2 





- r T i 


+ E 


J \Zl' x ' a \ 2 ds 







< oo . 



Such a solution exists under (HI), Lemma l2.1l and E[|C^ a: ' | 2 ] < oo. Note that taking the 
expectation in (|2.7p . Y^' x ' a = J(t,x,i;a). 
We now define for K > 0, 



a. (z A 



t.i 



E 



sup | Y, 

s£[t,T] 



t,x,a I 



< k(i + \x\ 2 )], 



and claim that for some constant K, the supremum in Vi(t,x) may be taken over a E 



A^'i (x). First taking the conditional expectation in (j2.T|) . we have 

< v h {X*?' a ,lZ) < K(l + IX^I), t < u < T, 

so that by Lemma l2. li the only restriction is to have a lower bound on Yu' x ' a . As in Lemma 
12.21 this is done by considering strategies with fewer interventions. Given a 6 <A?'J X) , consider 
the stopping time 

r = inf{s > t : J(s,X*' x ' a ,/ s Q ;a°) > Y s *' x ' a } 
where a is the strategy with no switches, and define a = (f n , t n ), where 

f n = T n l {rn < r} + OOl {Tn>r} . 

Now for each t < u < T, taking the conditional expectation in (I2.7P we obtain 

rt.x.a \H,x,a\ 



1 , , IV 1 ^^ _ V'^,«^ 

E 1 



'1' 



E 



fix**", I s )ds - g(X^' a ,I T ) + C*?' a - Ci' x ' a ) | T u 
l { n<.<T}(J(r,^,/ r a ;a ) - ^'«) | T u 



where we have taken the conditional expectation w.r.t. JF T inside the expectation. Since 
the process ( J(u, Xu X ' a , I"; a ) — Y u ' x,a \ < <T has right-continuous paths , by definition of 

r we have J(r, X* T ' x ' a , J? , a ) - Y" T M '' Q > a.s., so that 



1 {U < T} {Y^ & -Yt x n > 
Noting that on {u < r} we have 



(2.t 



it 1/ — "I - ^-**-n 



> J{u,X^ ,I u _;a u ) +c(X l ; ,I U -,I U ) 

> -K(1 + \X U \), 
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and since on {u > r}, Yu' x ' a = J(u, Xu X ' a , /„; a ), from Lemma 12.11 it follows that a E 



rt,x,a ja. n 



A s t ^(x), for some if not depending on (t,x). Furthermore taking u = t in (|2.8p . we have 
J(t,x,i;a) > J(t,x,i;a), and this proves the required assertion. 

• Step 4- Finally we show that for each K, there exists some positive K s.t. Al\ (x) C 
Agi(x). We fix a G ^'f (x). Applying Ito's formula to ly*^ 2 in ([22]), we have 



|y/' x '°| 2 + / |2£ x - a | 2 da = \g(X^ x ' a ,I^)\ 2 + 2 Yt' x > a f(Xl> x > a ,I°)d 
Jt Jt 



2 I Y*' x > a Zt> x ' a dW, - 2 I Y^dC 1 ^ 



Using (HI) and the inequality 2ab < a + b for a,i £ K, we get 

cT 

\Zl> x > a \ 2 ds < K{\+ sup \X l s ^ a \ 2 + 

se[t,T] se[t,T] 



f \Zt> x > a \ 2 ds < k(i+ sup |A^< Q | 2 + sup \Y*' x ' a \ 2 + \Cf£ x ' a - Ct' x ' a \ sup {Y^l) 



se[t,T] 

rT 

-2 I Y*' x ' a Zl' x ' a dWs 



(2.9) 



Moreover, from (j2.7|) . we have 



_ c t, X , a] 2 < + gup | X M,«|2 + gup |y*,*,a|2 



+ 



sG[t,T] 
T 



se[t,T] 



Zi' x ' a dW, 



(2.10) 



2 l2 

Combining (|2.9p and (|2.10p and using the inequality ab < ^ + for a, 6 G R and e > 0, 
we obtain 

/ |Z£ x ' a | 2 da < A-f(l + e)fl+ sup |X*' X ' Q | 2 ) + sup |Y s *' x ' a | 2 (e + -) 



+ e 



se[t,T] 

J Z^ a dW s 2 ) - 2 J Y*' x ' a Zl> x > a dW s . 



Taking the expectation in the previous estimate, it follows from Lemma \2. II and a 6 A s t '^ {x) 
that 



E 



|Z*' x ' a | 2 ds 



< K((l + e)(l + M sup |X* AQ | 2 
v v se[t,T] 



+ (e+-)E sup \Y. 

£ se[t,T] 



t,x,a\2 



+ eE 



Zt' x ' a dW, 



< K^{l + \x\ 2 )(l + e + -) +eE 



r 



^' x ' a | 2 (is 



Taking e small enough, this yields 

e' 



\Zl> x > a \ 2 ds 



'-Jt 



< K(l + \x\ 2 ), 
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Taking the expectation in (|2.10p . and using the previous inequality together with Lemma 
12.11 and a £ A^.\ (x), we get: 

E |^*,a* _£*,*,<*• |2 < K( 1 + | x |2) ; (2.11) 

for some positive constant K not depending on (t,x,i). Since (r n ) is strictly increasing, 
we know that at the initial time t, there is at most one decision time t\. Thus, from the 
linear growth condition on the switching cost, E[|C*' X ' Q | 2 ] < K(l + \x\ 2 ), which implies with 
(|2XTj) that a e A?i(x), and this proves the required result. □ 

In the sequel of this paper, we shall assume that (HI) and (He) stand in force. 

3 Time discretization 

We first consider a time discretization of [0, T] with time step h = T/m < 1, and partition 
= {tk = kh, k = 0, . . . , m}. For (t k , i) £T/,x I q , we denote by A% i the set of admissible 
switching controls a = (r n , i n ) n in At k) %, such that r n are valued in {£h, t = k, . . . , m}, and 
we consider the value functions for the discretized optimal switching problem: 

m— 1 

v?(t k ,x) = sup ¥.\Y J f{Xl^ a J tl )h + g{Xl^ a J tm ) 

JV(a) 

- ^c^'^in-i,^)], (3.1) 

n=l 

for € T h x I q x M d . 

The next result provides an error analysis between the continuous-time optimal switch- 
ing problem and its discrete-time version. 

Theorem 3.1 There exists a positive constant K (not depending on h) such that 

\vi{t k) x)-v?{t k ,x)\ < K(l + \x\ 5 / 2 )(hlog(2T/h)) 1/2 , 

for all (t k ,x,i) E^xR^x I q . 

If the cost functions Cij, i,j G I q , do not depend on x, then 

\vi(t k ,x)-vHtk,x)\ < K{1 + \xf 2 )h l f 2 



Remark 3.1 For optimal stopping problems, it is known that the approximation by the 
discrete-time version gives an error of order h% , see e.g. [12] and pQ. We recover this rate 
of convergence for multiple switching problems when the switching costs do not depend on 
the state process. However, in the general case, the error is of order (/ilog(l//i)) 2 . A rate 
of h2~ £ was obtained in [5] in the case of uncontrolled state process X, and is improved 
and extended here when X may be influenced through its drift and diffusion coefficient by 
the switching control. 



11 



Before proving this Theorem, we need the three following lemmata. The first two deal 
with the regularity in time of the controlled diffusion uniformly in the control, and the third 
one deals with the regularity of the controlled diffusion with respect to the control. 



Lemma 3.1 There exists a constant K such that 



sup max 



sup \X$™*-X%™*\ 
se[t e ,t£ +1 ] 



< K(l + \x\)h 7 - 



for all x G R d , i G L, k = 0, . . . , n. 



Proof. From the definition of X l ' x ' a in ([13]), we have for all (t k ,x,i) G T h x R d x l q and 

a G A tk ,i, 

E 



I yt,x,a yt,x,a\2 

sup |A U ' - A te 

ue[t e ,s] 



< K E 



\b Iu (X^)\du 



+ E 



sup 

u£[ti,s 



a Ir (X^ a )dW, 



for all s G From BDG and Jensen inequalities, we then have 



E 



sup \A U ~ A t( 

L u£[t e ,s] 



< k(e 



+ E 



t, 



\a Iu (X^)\ 2 du]) , 



From the linear growth conditions on 6j and Uj, for i £ I q , and Lemma 12.11 we conclude 
that 



E 



sup \xl> x > a - X*' x ' a|P 



s£[t e ,t l+1 ] 



< K p (l + \x\ p )h. 



□ 



Lemma 3.2 There exists some positive constant K such that 



sup 



sup \X lk 

0<s,u<T 
\s— u\<h 



t k ,x,a j£tk,x,a 



< 



K(l + \x\)(hlog(2T/h)) 



Proof. This follows from Theorem 1 in [7], using the estimates from Lemma 12. II and linear 
growth of bi, (7j. □ 

For a strategy a = (r n ,,t n ) n G -At fc ,i we denote by a = {f n ,l n ) n the strategy of A% i 
defined by 

f n = min{t£ £T t : t t > r n } , Z n = t n , n G N. 

The strategy d can be seen as the approximation of the strategy a by an element of A\ k i . 
We then have the following regularity result of the diffusion in the control a. 

Lemma 3.3 There exists some positive constant K such that 



sup \X Zk 
se[t k ,T] 



tl.,X,Ct \rtk,X,Ct 



1 

< K(E[N(a) 2 ]Y{l + \x\)h^ , 



for all x G M. d , i G l q , k = 0, . . . , n and a G A 



tk,i- 
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Proof. From the definition of X l ' x ' a and X*'*' 01 , for (t k ,x,i) G T h x R d x I q , a G .4^ 



we have by BDG inequality: 
E 



sup \X^ x ' a - X^ x > & \ 2 

u£[t k ,s] 



+ E 



\b(X^ a ,I u )-b(X^ & J u )\ 2 du 



tk 



\a{X^ a ,I u )-a{X^ & X)\ 2 du 



for all s E [tfc,T]. Then using Lipschitz property of 6j and <jj for i £ I q we get: 

1 2" 



E 



sup \Xl^ a - Xf> x 

u£[t k ,s 



< K[E 

+ E 
+ E 

< k(e 

+ E 



th 



tk 



tk 



\xt' x ' a - xh x ' & \ 2 du 



\b(X^ a ,I u )-b(X^ a ,I u )\ 2 du 



\a(X^ a ,I u )-a(X^ a ,I u )\ 2 du 



sup \X^ a - X^' 2 



da 



ue[t k ,T] 

for all s G [ifc,T]. From the definition of a we have 



tk re[tk,v] 

sup |X^| 2 + 1) / l Is , Ie ds 



tk 



•'tk 

which gives with (j3.2|) . Lemma 12.11 Remark 12.11 and Holder inequality: 



E 



sun I yt,x,a yt,a 
SUp |A U - A u 

«e[tfc,s] 



2" 


< If 


/ sup | 






-Jt k r£[t k ,u] 



rt,x,a 



>-t,x,a 1 2 



+ (E[N(a) 2 ])*(l + \x[ z )h), 



(3.2) 



□ 



for all s S [ifc,T]. We conclude with Gronwall's Lemma. 

We are now ready to prove the convergence result for the time discretization of the optimal 
switching problem. 

Proof of Theorem 13.11 We introduce the auxiliary function v\ defined by 



sup E 



N(a) 



f(Xi** a J 8 )d8 + g(X%>*> a ,I T ) - c(X% 



** n=l 

for all (t k ,x) £T h x R d . We then write 

\vi(t k ,x) -v?(t k ,x)\ < \vi(t k ,x) -v^(t k ,x)\ + \$(t k ,x) -vf(t k ,x)\ 
and study each of the two terms in the right-hand side. 
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• Let us investigate the first term. By definition of the approximating strategy a = (f n , Z n )n 
G J^l j of a E -At k> ii we see that the auxiliary value function v\ may be written as 



Vi(tf~,x) = sup E 



N(a) 



f{X$*> & j a )d8 + g{X%* & ,i T ) - c(4*'*'VVi, 



n=l 



where I is the indicator of the regime value associated to a. Fix now a positive number K 
s.t. relation (|2.6p in Proposition 12.11 holds . and observe that 



sup E 



N(a) 



f{Xy^X)ds + g{X^ & J T ) - J2 c(X^ & ,t^ u \ 
< Vi(t k ,x) < Vi(t k ,x) 



tk 



n=l 



sup E 



N(a) 



tk 



f(Xl^ a , I s )ds + g(X^ a , I T ) - £ c(X^> Q , t n ^,L n ) 



n=l 



We then have 



\vi(t k ,x) -v%(t k ,x)\ < sup Al (a) + A? (a) 



(3.3) 



with 
l 



E 



E 



If 




N(a) 

[E 

n=l 





Under (HI), and by definition of a, there exists some positive constant K s.t. 



A] k Ja) < K{ sup E 

< k( sup E 

y se[t k ,T] 



X 



+ E 



sup \X^\+l) / l^ds 

'se[t k ,T] 



+ 1 + 



sup |X* fe ' x ' Q 



( 



E 



L Jt 



by Cauchy-Schwarz inequality. For a 6 ^4f; j(x), we have by Remark |2. II 



E 



< KE 



N(a) < rjK^l + \x\)h, 



for some positive constant r] > 0. By using this last estimate together with Lemmata 12.11 
and 13.31 into (|3,4p . we obtain the existence of some constant K s.t. 



sup 



A* (a) < K(l + |x| 3 / 2 )^, 



(3.5) 
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for all (t k ,x,i) G T h x R d x I g . 

We now turn to the term A 2 x (a). Under (HI), and by definition of a, there exists some 
positive constant K s.t. 



N(a) 

aL» < *e[£ \ x % x,a - x t' xA] 



< K E 



n=l 
Ma) 



+ E 



n=l 

AT(Q!) 

'e[ |x! fe - x '' a _X>' X ' Q | 



N(a) sup IX^-X^ll) 



n=l 



iV(a) || 


sup X* fc ' x '° -Xl k ' x ' & \ 




s£[t k ,T] 



(3.6) 



by Cauchy-Schwarz inequality. For a G .4.^ j(x) with Remark 12. H and from Lemma 13.31 
we get the existence of some positive constant K s.t. 



N(a) 




sup 




2 


se[t k ,T\ 



I y\. s yv s 



< K(l + \x\ 5/2 )h^ . 



(3.7) 



On the other hand, 

N(a) 



E[El^ Q -4r 

n=l 



< E 



iV(a) sup \xl^ a - X^' x ' a \ 

0<s,u<T 
\s—u\<h 



< 



N(a) 



sup 

0<s,u<T 
|s— u\ <h 



\^s 



by Cauchy-Schwarz inequality. For a G .A^^x), by Lemma 13.21 this yields the existence 
of some positive constant K s.t. 

_N(a) 



n=l 



< 



if(l + |x| 2 )(/ilog(2T//i)) 1/2 



(3.8) 



By plugging ()3.7p and (|3.8p into (|3.6p . we then get 

A 2 » < A'(l + |x| 2 )(/ i log(2T//z)) 1 / 2 . (3.9) 
Combining (|3.5p and (|3.9p . we obtain with (|3.3p 

\vi{t k ,x)-v?(t k ,x)\ < K(l + \x\ 2 )(hlog(2T/h)) 1/2 . 
In the case where c does not depend on the variable x, we have A 2 x (a) = 0, and so by 



|ui(t fe ,x) -t£(t fc ,x)| < A'(l + |x| 3/2 )/ii 
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For the second term, we have by definition of and v^: 



m-i , tt+1 

%(t k ,x)-v?(t k ,x)\ < sup E ^ fiXl'-\L) f(X;f".l. 

^t,i i=k Jt * 



ds 



since I s = It e on [t^,t^ + i). Under (HI), we get 

\v!i (t k , x) — v^(t k , x)\ < K sup max sup E 

aeA u k<i<m-l se[te ,t e+1 ] 



for some positive constant K, and by Lemma 13. H this shows that 

\v?(t k ,x)-i#{t k ,x)\ < K{l + \x\)h*. 



□ 



In a second step, we approximate the continuous-time (controlled) diffusion by a discrete- 
time (controlled) Markov chain following an Euler type scheme. For any (t k ,x,i) £T/,x 



I q , a £ At j, we introduce {X^ tk,x,a )k<i< m defined by: 



X 



tk 



X, I''" 



h,t h ,x,a _ rph ( v h,t k ,x,a 



k < £ < m - 1, 



where 



F^x^k+i) = x + bi(x)h + Vi(x)vh fi k+1 , 



and "&k+\ = (Wt k+1 — Wt k ) / y/h, k = 0, . . . , m — 1, are iid, A/"(0, /^-distributed, independent 
of J~t k . Similarly as in Lemma \2.\\ we have the L p -estimate: 



sup 



max X 

£=k,...,m 



h,t k ,x,a I 
t, 



< KJl + \x\), 



(3.10) 



for some positive constant K p , not depending on (h,t k ,x,i). Moreover, one can also derive 
the standard estimate for the Euler scheme, as e.g. in section 10.2 of jllj : 



sup 



I Y l k,x,a Y h,t k ,x,a\ 

maX A t< ~ ^t t 
=k,...,m 



< K p (l + \x\)Vh. (3.11) 



We then associate to the Euler controlled Markov chain, the value functions v^, i € I q , for 
the optimal switching problem: 



m— 1 



Vi(t k ,x) 



sup E\Y,f(X?f k ' X ' a Jtt)h + 9(xl 
l=k 



tu ,x,a 



It m ) 



N(a) 

^ c (X^Vn-l^n) 



(3.12) 



n=l 



The next result provides the error analysis between v\ by vf, and thus of the continuous 
time optimal switching problem Vi by its Euler discrete-time approximation v^. 
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Theorem 3.2 There exists a constant K (not depending on h) such that 

\v?(t k ,x)-v?(t k ,x)\ < K{l + \x\ 2 )Vh, (3.13) 
for all (t k ,x,i) G T h x M d x I q . 

Remark 3.2 The above theorem combined with Theorem 13.11 gives the rate of conver- 
gence for the approximation of the continuous time optimal switching problem by its Euler 
discrete-time version: there exists a positive constant K s.t. 



\ Vi (t k ,x)-4(t k ,x)\ < K(l + \x\ 5 / 2 ){hlog(2T/h))K 



(3.14) 



for all (t k , x, i) 6 Tj x R** x I q . Moreover if the cost functions Qj, i,j G I q , do not depend 
on x, then 

\vi(t k ,x)-v?(t k ,x)\ < K(l + \x\ 2 )h^, 
Proof of Theorem 13.21 

• Step 1. For (t k ,x,i) G T/ 4 x R d x I q , and a G A\ i we denote by 



m— 1 



N(a) 



J h (t k ,x,i;a) = Efj/^^ + ^IJ- E c (^' a > 



'w— 1 1 'r, 



71=1 



so that Vi(t, k , x) = sup^g^h J h (t k ,x, i, a). Given a G i , let us define i 7 " = f(X t *' x ' a , I" e ), 

c a = c^^J^,/") and = E[ 1)^(^-^)1^], for I = k,...,m. Consider 
the stopping time 

r = M{t t >t k : J h (t e ,X%> x > a ,I? t ;a°)>Y?}, 
where a is the strategy with no switches, and define a = (f n , t n ), with 

f n = T n l {Tn < r} + ool {Tn>r} . 



As in the proof of Proposition 2.1, we easily check that 
and 

Y? > J(tt,X%**I*;aP), 



(3.15) 



(3.16) 



for all £ = k, . . . , m. From (|3.16p and the estimates on X^' x ' a in Lemma [2. H we know that 



E 



sup (|Ff| 2 + \F?\ 2 + \cf\ 2 ) < K{l + \x\ 2 ), 

k<l<m - 1 



(3.17) 



for some positive constant K. Moreover, by definition, we have: 

if = E [Y e & +1 \T te ] +hF e -c e , £ = *,..., ro - 1. 
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Letting AM" +1 := Y^_ 1 — E[Y^ 1 |J r t J, we obtain in particular 



m— 1 



m— 1 



m— 1 



5>f = ft^Ff-^AM^ + ^-lf), 



i=k 



and so by f|3. 17j) 



E \^24\ < A'(l + |x| 2 ) + 3E 
e=k 



'm—1 



E AM e+i 

i=k 



m—1 



i=k 



m—1 



= K(l + \x\ 2 ) + 3E 

Now by writing that 

m—1 



(3.18) 



£=fc 
m—1 



E (AMf +1 - + cf )(23f + AMf +1 - hFf + c 



t ) > 



we get 



m—1 



m—1 



m—1 



1=0 



^|AM^| 2 = |y«| 2 -|y «| 2 -^/ l Ff(/ l F/-2F/-2cf)-2^ C fy, 

^=fc £=0 

m—1 m—1 

- ^Mf +l {2Y? - 2hFl + 2cf) - E kf 
e=o 

Since E [AMf +1 | = 0, this shows 

m—1 

~\Y£\ 2 - hF?{hF? - 2Y? - 2cf) - 2 E 4rf 



t=o 



E[E|AMf +1 p 



m—1 



m—1 



e=k 



< E 



1=0 



e=o 



m—1 



< K{1 + |x| 2 ) + 2E [| E c f y / 



(3.19) 



e=o 



where we used again (|3.17p . Now since ct > 0, 



E 



m—i m—l 

|E c f y /l ^ 4(E c ^) su p i y ^ 5 

Ll ^ L V ' / k<£<m-l 



e=o 



e=o 

m—1 

< eE[^|AM £ % 



+ K(1 + -)(1 + 



X 



for all e > 0, by (|3.17p . (|3.18p and Cauchy-Schwarz inequality. Hence taking e small enough 
and plugging this estimate into (|3.19p . we obtain 



m—1 



e[^|AM,%| 2 ] < K(l + \x\ 2 ). 
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Using (|3.18p one more time and recalling that N(a) < r/ *£2e cf for some r] > under the 
uniformly lower bound condition in (He), we thus obtain 

E|iV(a)| 2 < A'(l + |x| 2 ). (3.20) 

Combining this last inequality with (|3.15p , we get that the supremum in the definition (|3.ip 
of v£(t k ,x) can be taken over A^(x) = {a £ A% k>i s.t. E\N(a)\ 2 < K(l + |x| 2 )}. Using 
the same argument with X tk,x,a instead of X tk,x,a and estimate (|3.1U|) on ||X t ' ! '' tfe ' :r ' a || we 



also get that the supremum in the definition ([37L2]) v^(t k ,x) can be taken over A t ' i (x). 
• Step 2. Now, for any a G A^ i (x), we have under (HI) and by Cauchy-Schwarz inequality 

m— 1 



E 



[ £ h\f(X^ a J te ) - f(X^ a ,I te )\ + \g(Xt^ a ,I tm ) - g(X^ a ,I tm )\ 
t=k 
N(a) 

+ £ \c(X^ a ,L n ^,L n )- C (X^ a ,L n ^,^ 
n=l 



< K¥\(l + N{a))( sup \xl*' x ' a - X% tk ' x ' a \) 

k<e<m 1 



< K(l + \x\ 



sup \X t 



X 



k<l<m 

< K(l + \x\ 2 )^h, (3.21) 
by (|3.1ip . Taking the supremum over a € A^^(x) into (j3.21|) . this shows that 

\v?(t k ,x)-v?(t k ,x)\ < K(l + \x\ 2 )Vh. 



□ 



4 Approximation schemes by optimal quantization 

In this section, for a fixed time discretization step h, we focus on a computational appro- 
ximation for the value functions v^, i € I q , defined in (j3. 12|) . To alleviate notations, we 
shall often omit the dependence on h in the superscripts, and write e.g. vi = v^. The 
corresponding dynamic programming relation for Vi is written in the backward induction: 

Vi(t m ,x) = gi(x), 

Vi(t k ,x) = max ^E[vi(t k+1 , Xj^ 1 )] + fi(x)h , max[vj(t k , x) - %(a;)]}, 
for k = 0, . . . , m — 1, (i, x) € I q x M. d , where X tk,x ' 1 is the solution to the Euler scheme: 

= F h^k+i) ■■= x+biWh+ai^Vhtik+L 

Observe that under the triangular condition on the switching costs Cjj in (He), these 
backward relations can be written as an explicit discrete-time scheme. Indeed, if Vi(t k ,x) = 
Vj(t k ,x) — Cij{x) for some j / i, we then have for / ^ i,j, 

Vj(t k ,x) - Cij(x) = Vi(t k ,x) > vi(t k ,x) - cu(x) 

> Vl(t k ,x) - Cij(x) - Cjl(x), 
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so that Vj(tk,x) > vi(tk,x) — Cjiix). By positivity of the switching costs, we also have 

Vj (t k ,x) =Vi (t k ,x) + Cij(x) > Vi (t k , X) - Cji (x) . 

It follows that 

Vj(t k ,x) = E[^(t fc+1 ,X^)] +fj(x)h, 

and (recalling that ca(-) = 0), the backward induction may be rewritten as 

Vi(t m ,x) = gi(x) (4.1) 
Vi(t k ,x) = max \E[vj(t k+1 , X^ 3 )] + fj^h-Cijix)}, (4.2) 

for k = 0, . . . , m — 1, (i,x) G I q x M, d . Next, the practical implementation for this scheme 
requires a computational approximation of the expectations arising in the above dynamic 
programming formulae, and a space discretization for the state process X valued in 
We shall propose two numerical approximations schemes by optimal quantization methods, 
the second one in the particular case where the state process X is not controlled by the 
switching control. 



4.1 A Markovian quantization method 

Let X be a bounded lattice grid on M. d with step 5/d and size R, namely X = {5/d)Z d n 
B(0,R) = {x G R d : x = (S/d)z for some z G Z d , and \x\ < E}. We then denote by Proj x 
the projection on the grid X according to the closest neighbour rule, which satisfies 

|s-Proj x (x)| < max(|s| - R, 0) + <5, Vx G R d . (4.3) 



At each time step t k G T^, and point space-grid x G X, we have to compute in (|4.2[) ex- 
pectations in the form E (p(X^f^*) , for </?(.) = v^{tk+i, i G I q . We shall then use an 
optimal quantization for the Gaussian random variable i?fc+i, which consists in approxima- 
ting the distribution of i9 ~» M(0, Id) by the discrete law of a random variable i? of sup- 
port N points wi, I = 1,...,N, in M. d , and defined as the projection of $ on the grid 
{wi, . . . ,wn} following the closest neighbor rule. The grid {w\, . . . ,wn} is optimized 
in order to minimize the distorsion error, i.e. the quadratic L 2 -norm — i?|| 2 . This 
optimal grid and the associated weights {tti, . . . ,ttn} are downloaded from the website: 



http:/ /www■quantize■maths-fi■Com/downloads| , . We refer to the survey article [15] for 



more details on the theoretical and computational aspects of optimal quantization meth- 
ods. In the vein of |16j . we introduce the quantized Euler scheme: 

X%> + f = Proj x (ifM)), 

and define the value functions Vi on T m x X, i G I q in backward induction by 

Vi(t m ,x) = gi(x) 

Vi(t k ,x) = max lE[vj(t k+1 , X^ff)] + fj(x)h -Cij(x)\, k = 0, . . . , m - 1. 
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This numerical scheme can be computed explicitly according to the following recursive 
algorithm: 

Vi(t m ,x) = gi(x), (x,i)£XxI q 
N 

Vi{t k ,x) = max y^7T/ v j (t k+ i,Pioi x (Fj l (x,wi))) + fj(x)h - Cij(x) , (x,i) GXx I q , 

J& g L j=i 

for fe = 0, . . . , m — 1. At each time step, we need to make O(N) computations for each point 
of the grid X. Therefore, the global complexity of the algorithm is of order 0(mN(R/S) d ). 

The main result of this paragraph is to provide an error analysis and rate of convergence 
for the approximation of V{ by ■Oj. 

Theorem 4.1 There exists a constant K (not depending on h) such that 
\vi(t k ,x) -Vi(t k ,x)\ < Kexp(Kh~ 1 \\$-m)(l + \x\ + -^ 

^+"- 1/2 ll»-*(i+M + f 

+ J_ exp(A:ft -2||j_^ ) ( 1 + N 2 + ( ^ 

for all (tk,x,i) G T/, x X x I,. In the case where the switching costs Cij do not depend on 
x, the above estimation is stengthened into: 

\vi(t k ,x)-Vi(t k ,x)\ < K\h- 1 / 2 \\tf-$\\ 2 exp(Kh- 1 \\$-m)(l + \x\ + 



h 

Remark 4.1 The estimation in Theorem 14. II consists of error terms related to 



• the space discretization parameters 5, R, which have to be chosen s.t. S/h and 1/Rh 
go to zero. 

• the quantization error — of the normal distribution JV(0, Id), which converges 

to zero at a rate A^3, where N is the number of grid points chosen s.t. h~ 'N~ goes 
to zero. 

By combining with the discrete-time approximation error (|3.14p . and by choosing grid 
parameters 5, 1/R of order bA , and a number of points iV of order l/h d , we see that the 
error estimate between the value function of the continuous-time optimal switching problem 
and its approximation by Markovian quantization is of order h?. With these values of the 
parameters, we then see that the complexity of this Markovian quantization algorithm is 
of order 0{l/h 4d+1 ). 

Let us now focus on the proof of Theorem 14.11 First, notice from the dynamic pro- 
gramming principle that the value functions Vi, i G I q , admit the Markov control problem 
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representation: 



m— 1 



Vi(t k ,x) 



sup E\J2 KXlf'^ !u)h + 9{Xt x ' a , I 
'^t k ,i i=k 

N(a) 

-^c(X^ )in _ 1)in ) 



where X tk,x,a is defined by 



X, 



te+i 



n=l 



Proj x (fJ (1J'*'V, +1 )), fc < € < m- 1, 



(4.4) 



for q 6 .4^ i; and A; = 0, ... ,m — 1, are iid, ^-distributed, and independent of T\ 



tk ■ 



We first prove several estimates on X tk ' x ' a . 

Lemma 4.1 For each p > 1 there exists a constant K p (not depending on h) such that 



sup 

i&A^ ,,fc<^<m 

z k> 1 ~ ~ 



X 



tk > x >& 



+ 



< 



sup 



(4.5) 



K p exp (K p h- p / 2 \\#-3\fy (l + \x\ + ^\ , 



for all (tk,x,i) £ T^xXxI,. 



Proof. We fix (t k ,x,i) eT h xXxI q , ae A% k>i , and denote X tl = X^' x ' a , k < i < m. 
Denoting by E; the conditional expectation w.r.t. J 7 ^, by a standard use of Gronwall's 
lemma and linear growth of bi, o~i, we have 



< e Kph 



X 



+ K v h. 



We will use the following convexity inequality : for a, b £ R+, h G [0, 1], 

(a + hbf < (1 + K p h)a p + K p hiP . 
By definition of F h , and the fact that |Proj x (?/)| < \y\ + 5 for all y £ R a 



(4.6) 



(4.7) 



X, 



te+i 



< 



FUx h ^ l+1 ) + h 1 ^a Jt (X tl )\S e+1 - t+1 \ + 5 



Fj h H (X t£ ,^ +1 ] 



+ h 



h 1 / 2 



+ 



h 



Combining this last inequality with (|4.6p . (|4.7p . linear growth of cij and the fact that 
are independent of Tt t -, we obtain 



< (1 + K p h) {e Kph \X te \ p + K p h) + A> 



' ^(x f ,)||^-^||; ^ 

hvl 2 + ftp 



1 + p _ #\\ p p h- p / 2 + ^ 
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By induction, taking the expectation, recalling that h = and since (l + ^) m < e y for 
all y > 0, we obtain 



E 



< /^exp^-f/ 2 ^-^^) (^l + \ x \P + ^- + h-v/ 2 p-$\\ p ^ 



for all k < £ < m. The estimate for F h (Xt l , , d£ + i) then follows from (|4.6p . 
Lemma 4.2 There exists some constant K (not depending on h) such that 



□ 



sup 



'ft* 1 



aim I y'*" 1 ' Y tk,x ' a \ 
SU P ~ A io 



k<e<m 



< K 



h- l/2 U-$\ 



1 



,exp 



(Kh- 1/2 p - S\\ 2 )(l + \x\ + 



+ 71 + Rh exp II* " <€) I 1 + M + (*) 



(4.8) 



/or aZZ (tk,x,i) G T/j x X x I g . 



For fc < ^ < m — 1, we get 



Proof. As before we fix (ti.,x,i), a and omit the dependence on (tk,x,i,a) in Xt e . Let us 
first show an estimate on 

On the other hand, since 



X te+1 -F? te (X te ,V e+1 ) 2 + F? tt (X tt ,0 i+1 )-Ff u (X tt ,-& t+1 ) 



F? t (X ti ,0 t+ i)-F? t (X tt ,# e+1 ) 



2 

(4.9) 



\v\ 2 

|y-Proj x (y)| < 5 + \y\l M > R} < S + 



by inequality (|4.3p . we have 



< 5 + 



<(X £ ,^ +1 ) 



(4.10) 



2 

Furthermore by standard estimates for the Euler scheme (see e.g. Lemma A.l in [15]), we 
have 



F?(X tt ,# t+1 )-Fk(X tt ,# i+1 ) n < (1 + Kh) 



Xt„ — Xt f 



and by the linear growth property of a and the fact that r &£ + i,'&£ + i are independent of J r t e , 



FUx tl ,di +1 )-FUx tv # l+1 ) 



< Kh 1/2 ( 1 + 



Plugging these three inequalities into (|4.9p . we get : 



X te+1 - X te+1 



< (1 + Kh) 



X te - X tl 



Kh 1 ' 2 



Xi 



x± 



|tf-tf|L (4.11) 



+ 1 



+ 5 + 



FUx t J i+1 ) 



R 
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Finally since X tk = X tk = x, we obtain by induction, and using the estimates f|4.5j) on 
Fl(X t J i+1 ) 



Xt„ — Xt 



< K 



/i- 1/2 ||^-i?|| 2 exp {Kh- x \\$-$\Q[l + \x\ + -) + - 



6 



Rh 



exp (Kh~ 2 p - (l + \x\ 2 + (-) 



(4.12) 



for all k < £ < m. Now by definition of Xt k , Xt k , we may write for k < £ < m — 1: 

Xt t+1 — X te+1 = (X te — X tt ) + h[b(X te , It e ) — b(X te , It e )) 
+ Vh{a(X te ,I te )$ £+1 - *(X ti ,I ti )0t+i) 
+ Projx « {X tl , &+l)) - Ft (X te , m ) , 



Since Xf k = Xt k (= x), we obtain by induction: 

m— 1 



sup 

k<Km 



Xt, — Xt 



e=k 

+Vh\\ sup | y^jr(X tr , it r )flr+l ~ Q'(X tr ,I tr )<& r+ i\ 

k<£<m r<£ 

m—1 

+ £|Proj x (fJ (X ti A + i)) ~ Fl{X t J t+l ) (4.13) 



l=k 



We now bound each of the three terms in the right hand side of (|4.13p . First, by the 
Lipschitz property of b and ()4.12j) . we have 

m— 1 

hJ2\Hx ti ,it,)-Kx ti ,i ti )\\ 2 



£=k 



< K\h~ 1/2 p - t?|| 2 exp {Kh- X p - $f 2 ) (l + \x\ + - 



1 



+ - + —exp(Kh- 2 \\$-m)(l + \x\ 2 + 



h J 



Next, recalling that is independent of J r t e , with distribution law d, and since # is an 
optimal L 2 -quantizer of it follows that E[^ +1 | F h ] = E[i9] = E[i?] = 0. Thus, the process 
(52r<£ a (^trihr)^r+i - o{X tr ,I tr )'& r j r i)e, is a 7"^ -martingale, and from Doob's inequality, 
we have: 



SUp | yVpQ^itJffr+l - <7(^tr ,it r )^r+l 

:<i<m r<£ 

m—1 



H+l\ 



By writing from the Lipschitz condition on <7j that 



\a(X te , / t< )0* +1 - a(X t£ ,I tt )0i+i\ < K [\ X U ~ X U \\ &t+i 

+ (1+|^| 2 )|^ 



l+l — V£ + l 
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and since are independent of Tt t -, we then obtain 

Vh SUP \^2a{Xtr,Itr)^r+l-<T(Xtr,Itr)^ 



k<£<m. 



r+l\ 



< K sup \\\X tt -X tt \l + {l + \\X tl \\ 2 )U-^\ 

k<e<m-l 1 



< K 



h~ l / 2 \\d - t?|| 2 exp (Kh- 1 p - m) (1 + \x\ + - 

+ 1 + ik exp ( Kh ~ 2 P - i 1 + n 2 + ,K2 



h' 



where we used the estimates (14. 5f) and (14. 12h . Finally the third term in (|4.13p is bounded 
as before by (|4TT0|) . □ 

Proof of Theorem 14.11 For (t/., x, i) £ x X x I q , we show as in the proof of Theorem 
13.21 that we can restrict to strategies a € A^ k i such that 



E|jV(a)| < K[l+ sup 



X. 



t k ,x,a 



for some constant K, not depending on (tk,x,i,h). By using the estimation (|4.5p . this 
means that the supremum in the representation (|3.ip of Vi(tk,x) can be taken over the 
subset 

j}fi(x) = {a£ A^iS.t. E|Af(a)| 2 < ^exp^/i- 1 )^-^)^) ( X + M 2 + }• 

Then, for a € ^4^^ (a;), we have under (HI) and by Cauchy-Schwarz inequality 
m— 1 

E[ £ fc|/(XJ'*>V*J - /(*£ ,a!, Vt,)| + \g(Xt^J tm ) - g(Xi^ a J tm )\ 

N(a) 

+ ]T ^{X^^^in)-^^,^^ 



n=l 



< KE 



[(l + JV(a))( sup Ix^'^-X^'^l 



k<£<m 



< K exp ( Kh~ l \\tf -tiWl) {l + \x\ + 



sup \X t - X t 

k<l<m 



< Kexp(Kh-^ti-^)(l + \x\ + --) ^ + h" l/2 p-^\\ 2 {l + \x\ + ^ 



(4.14) 



by Lemma l4.2i Taking the supremum over a G ^4^^"(x) in the above inequality, we obtain 
an estimate for \vi(tk,x) — Vi(tk,x)\ with an upper bound given by the r.h.s. of (j4. 14j) . 
which gives the required result. 
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Finally, notice that in the special case where the switching cost functions Cjj do not 
depend on x, we have 



m— 1 



\vi(t k ,x) - Vi(t k ,x)\ < sup E\^2h\f(xl 

~ Ah L 



k ,x,a 



a&A t k ,i £=k 



+ \g(&™I tm )-g(X%*~,I u 



< K sup E|X 



aeA*} ,k<£<m 



< A" 



/i" 1 / 2 !!^ - i?|| 2 exp {Kh~ l p - $\\ 2 ) (1 + \x\ + - 



by the estimate in Lemma 14,21 



□ 



4.2 Marginal quantization in the uncontrolled diffusion case 

In this paragraph, we consider the special case where the diffusion X is not controlled, i.e. 
bi = b, <Ji = a. The Euler scheme for X, denoted by X, is given by: 

Xq = Xq, X tk+1 = F h (X tk ,'&k+i) 

:= X tk +b(X tk )h + a(X tk )Vh^ k+1 , fc = 0,...,m-l, 

where $k+l = (Wt k+1 — Wt k )/Vh, k = 0, . . . ,m— 1, are iid, J\f(0, /^-distributed, independent 
of Tt k . Let us recall the well-known estimate: for any p > 1, there exists some K p s.t. 



\X tk \\ < K p (l + \\X \ 



(4.15) 



\p — i"- ii "up' 

Notice that the backward dynamic programming formulae (|4.ip - (|4.2p for V{ can be written 
in this case as: 



Vi(t m ,-) = 9i(-), i G lg 

Vi(t k ,.) 



max[P Vj(t k+1 , .) + hfj - cy]. 



Here P is the probability transition kernel of the Markov chain X, given by: 
P h y{x) = E[cp(X tk+1 )\X tk =x] = E[ip(F h (x, ■&))], 



(4.16) 



(4.17) 



where •& is M(0, /^-distributed. Let us next consider the family of discrete-time processes 



t h )k=0,...,m 



i £ Iq, defined by: 



Vi(t k ,X tk ), k = 0, . . . ,m, i € Iq. 



Remark 4.2 By the Markov property of the Euler scheme X w.r.t. (J-t k ) k , we see that 
(Yi k )k=o,...,mi i S Ig> satisfy the backward induction: 



9i(X t 



9i{X T ), i € Iq 



K = ^{EK J^+hfjiX^-c^X^)}, k = 0,...,m-l, 
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and is represented as 

Yi = ess sup E I" f(X tt , + g (X tm , I tm ) - Yl c (^« ' ^ > ^ 



m— 1 



JV(a) 



n=l 



On the other hand, the continuous-time optimal switching problem (|2.4p admits a repre- 
sentation in terms of the following reflected Backward Stochastic Differential Equations 
(BSDE): 



9i (X T ) + [ T f(X s )ds - [ T Z\dW s + K l T - K\, i G lg, < t < T, 

T 



Yi > max[Y/ - CidXt)] and / (Y£ — max[Y/ — Cij(X t )])dKi = 0. (4.18) 

Jo j& 

We know from [6], [10] or [9] that there exists a unique solution (Y, Z, A') = (Y 4 , A l )j G n 9 
solution to (|4.18p with Y G 5 2 (M' 3 ), the set of adapted continuous processes valued in IR 9 
s.t. E[sup 0<t<T \Yt\ 2 ] < oo, Z G 7W 2 (1R' ? ), the set of predictable processes valued in W 1 s.t. 
E[J Q T |Z t | 2 <ii] < oo, and K l G S 2 (R), Aq = 0, K l is nondecreasing. Moreover, we have 



Y l t = Vi (t,X t ), iel q , 



ess sup E 



rp N(a) 

f(X s ,I s )ds + g(X T ,I T ) - Y 4X Tn , 



n=l 



, < t < T. 



We propose now an optimal quantization method in the vein of pQ for optimal stopping 
problems, for a computational approximation of (Y2 )k=o,...,m- This is based on results 
about optimal quantization of each marginal distribution of the Markov chain (-Xt fe )o<fc<m- 
Let us recall the construction. For each time step k = 0, . . . , m, we are given a grid 
= {x\, . . . ,x k h } of N k points in W 1 , and we define the quantizer X k = Pro] k (X tk ) of X tk 
where Proj^ denotes a closest neighbour projection on Tf~- For N k being fixed, the grid 
is said to be L p -optimal if it minimizes the L p -quantization error: \\Xt k — Proj fc (At fc )|| p . 
Optimal grids are produced by a stochastic recursive algorithm, called Competitive 
Learning Vector Quantization (or also Kohonen Algorithm), and relying on Monte-Carlo 
simulations of X tk , k = 0, . . . , m. We refer to [15] for details about the CLVQ algorithm. 
We also compute the transition weights 



nx k 



+i 



X, 



4] 



F[(X tk+1 ,X tk ) £ Q,(T k+1 ) x QjT,)] 

p[x tk g d(T k )] 



where Ci(T k ) C {x G M. d : \x — x l k \ = min^r*, \x — y\}, 1 = 1,..., N k , is a Voronoi tesselation 
of Fj.. These weights can be computed either during the CLVQ phase, or by a regular 
Monte-Carlo simulation once the grids T k are settled. The associated discrete probability 
transition P k from X k to X k+ \, k = 0, . . . , m — 1, is given by: 



E 

v=i 



E[ip(X k+1 )\X k 



41 
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One then defines by backward induction the sequence of IR^-valued functions vt = (^)iei, 
computed explicitly on T^, k = 0, ... ,m, by the quantization tree algorithm: 



% = max [PkV 3 k+1 + hfj - c^-] , k = 0,...,m-l. 



(4.19) 



The discrete-time processes (Y£ )k=o,...,m, i £ ar e then approximated by the quantized 
processes {Y k )k=o,...,m-, i G Iq defined by 

= fe = o,...,m, i e v 



The rest of this section is devoted to the error analysis between Y % and Y % . The analysis 
follows arguments as in [2] for optimal stopping problems, but has to be slightly modified 
since the functions •) are not Lipschitz in general when the switching costs depend on 
x. Let us introduce the subset LLip(M. d ) of measurable functions 99 on M. d satisfying: 

\<p(x) - <p(y)\ < K(l + \x\ + \y\)\x-y\, Vx,yeM d , 
for some positive constant K, and denote by 

\<p(x) - <f(y)\ 



[<P] 



LLip 



sup 



x. 



,x^ y (1 + \x\ + \y\)\x-y[ 



Lemma 4.3 The functions Vi(tk, ■), k = 0, . . . , m, i G I q , lie in LLip(M. d ), and [vi(tk, -)] LLip 
is bounded by a constant not depending on (k,i,h). 



Proof. We set v l k = Vi(tk, ■)■ From the representation (13.12j) . we have 

m-l N(a) 

> E| £/(X^Vt> + 0(^'VtJ- 



sup E Vt> + 0(^'VtJ- ^c(X^, tn _ 1)tri 

'^ A t,i £=k n=l 



where X tk,x is the solution to the Euler scheme starting from x at time tf.. From (|4.15p . 
we see, as in the proof of Theorem I3.2( that in the above representation for v k (x), one 
can restrict the supremum to A^^(x) = {a G A% kii s.t. E|iV(a)| 2 < K(l + |x| 2 )} for some 
positive constant K not depending on (tk,x,i,h). Then, as in the proof of Theorem 14.1} 
we have for any x,y G M. d , and a G ' i (x) U A t '^(y), 



m—l 



E 



[Y,h\f(xZ*M-f(Xu ,v M\ + \9(xt x JtJ-9(x£\it m )\ 



l=k 



N(a) 



t,,,x 



» = 1 



< K(l+\\N(a)\\ 

k<£<m 

< K(l + \x\ + \y\)\x-y\, 



clln I y^^ ytk,y\ 
SU P \ X t e ~ X t e I 



2S 



by standard Lipschitz estimates on the Euler scheme. By taking the supremum over A^^(x) 
U A^' K i (y) in the above inequality, this shows that 



\v\(x) - v\(y)\ < K(l + \x\ + \y\)\x-y\, 



i.e. v l k S LLip(U. d ) with [v 



□ 



The next Lemma shows that the probability transition kernel of the Euler scheme 
preserves the growth linear Lipschitz property. 

Lemma 4.4 For any (p £ LLip(M. d ), the function P h ip also lies in LLip(M. d ), and there 
exists some constant K , not depending on h, such that 

[P h ^ iP < V3(i + o(h)M LLw , 

where 0(h) denotes any function s.t. 0(h) jh is bounded when h goes to zero. 
Proof. From (|4.17|) and Cauchy-Schwarz inequality, we have for any x, y G M. d : 
\P h y(x) - P h ip(y)\ 

,2\V2 



< (E\<p(F h (x,0))-<p(F h (y,#))\ 

< M LLip (E| (1 + \F h (x, 0)| + \F h (y, $)\) 2 \F h (x, 0) - F h (y, ) | 2 ) ^ 

< v^M^fEfa + l^^^^ + lF^y^)! 2 )!^^^)-^^^)! 2 ]) 1 , (4.20) 



where we used the relation (a+b+c) 2 < 3(a 2 +6 2 +c 2 ). Since i? has a symmetric distribution, 
we have 



(1 + \F h (x,$)\ 2 + \F h (y,$)\ 2 )\F h (x,$) -F h (y,$)\ 2 
i (1 + \F h (x, $)\ 2 + \F h (y, tf)| 2 ) \F h (x, t?) - F h (y, $)\ 2 
+ (1 + \F h (x, -tf)| 2 + \F h (y, -$)\ 2 ) \F h (x, i?) - F h (y, ■ 



A straightforward calculation gives 

i[(l + \F h (x,$)\ 2 + \F h (y,$)\ 2 )\F h (x,$) - F h (y,$)\ 2 

+ (1 + \F h (x, -tf)| 2 + \F h (y, -$)\ 2 )\F h (x, -0) - F h (y, -$)\ 2 

= (l + \x + hb(x)\ 2 + \y + hb(y)\ 2 + h\o(x)$\ 2 + h\a(y)d\ 2 ) \x - y + h(b(x) - b(y))\ 2 
+ h\(a(x) - a(y))$\ 2 (\x + hb(x)\ 2 + \y + hb(y)\ 2 ) 

+ AhUx + hb(x)\a(x){f) + (y + hb(y)\a(y)0)] (x - y + /i(6(x) - 6(y))|(tr(s) - <x(y))0) 
+ h 2 (\a(x)$\ 2 + |a(y)0| 2 )|(a(x) - a(y))V\ 2 . 

By Lipschitz continuity of b and <r, and the fact that E|i9| 4 < oo, we deduce that 



E 



(1 + \F h (x,$)\ 2 + \F h (y,$)\ 2 )\F h (x,$) - F h (y,$)\'< 



< (l + 0(h))(l + \x\ 2 + \y\ 2 )\x-y\ 2 . 
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Plugging this last inequality into (|4.20p shows the required result. 



□ 



We now pass to the main result of this section by providing some a priori estimates for 
\\Yt k — Y k \\ in terms of the quantization error \\X tk — X k \\. 

Theorem 4.2 There exists some positive constant K , not depending on h, such that 

m 

\Yl k -n\\ p < KY,(l + \\Xo\\r + \\X e \\ r )\\X te -X e \\ s , (4.21) 



max 

iei q 



£=k 



for any k = 0,... ,m, and (p,r,s) £ (l,oo) s.t. | = ~ + i. 

Proof. We set v\ = Vi(t k , .), and by misuse of notations, we also set Y k % = Y t l k = v k (X k ). 
From the recursive induction (|4.16p (resp. (|4.19p ) on v k (resp. v k ), and the trivial inequality 
| max 7 - aj — max,' a A < max, |o ? - — a A, we have for all i 6 L: 



\n-Yi 



= \vi(x tk )-vi(x k )\ 

< max | [P h v{ +l (X tk ) + hf 3 {X tk ) - c l3 (X tk )] - [P k v{ +l (X k ) + hf 3 {X k ) - c l3 {X k ) 

'J&q 



< max 

j&g 



P h 4+i(Xt k ) - P k v{ +1 {X k )\ + h\f 3 {X tk ) - f 3 (X k )\ + \c l3 {X tk ) - c t3 (X k ) 
hs-J 



< K X tk _ Xk \+ max \P h v{ +l {X tk ) - P k v{ +1 (X k 

J&q 

by the Lipschitz property of /» and Cj,-, and so 



max 


\ri \ri 
x k ~ x k 


< K 


Xt k - x k 


+ max 


i&q 




P 


P i&q 



p h 4 +1 (x tk ) - p k vi +1 (x k 



(4.22) 



Writing E^ for the conditional expectation w.r.t. X k , we have for any j G I g 

\P h 4 + i(Xt k )-Pkvi + i(X k )\ 
< \P h 4 +1 (X tk ) - P h v{ +l (X k )\ + \P h vUi{Xk) ~ E k [P h vl +1 (X tk )}\ 



+ \E k [P h vl +1 (X tk )]-P k vl +1 (X k )\ 



P h 4 + i(Xt k ) - P n vl +1 (X k )\ + \E k [P%l +1 (X k ) - P n vl +1 (X tk )} 



jh-i 



h-i 



yh-i 



+ \E k [Y k i +1 -Y k i +1 }\. 
Since E^ is a L p -contraction, we then obtain 



< 2 



P h 4 +1 (X tk )-P h 4 +1 (X k 



v l ~V~i 

x k +i I k +i 



< K{\ + 0{h)) (l + \X tk \ + \X k \)\X tk -X k \ 



+ 



Y L 



fc+i 



I k +i 



< A'(l + O(fc))(l + ||X || r + ||**|| r ) 



X, 



x k 



+ 



* k +i ~~ * k +i 



(4.23) 



30 



where we used Lemmata 14.41 and 14.31 Holder's inequality and (|4.15j) . Substituting (|4.23|) 
into (f4T22|) . we get 



max 



I u — I u 



< K(l + 0(h))[l + \\Xo\\ r + \\X k \ 



Xu. — Xk 



+ max 

tela 



*k+l ~ x k+l 



for all k = 0, . . . ,m — 1. Since maxj 6 i II — Y£ 



max ieI \\gi(X tm ) - g(X n 



K\\Xtm — X m \\ by the Lipschitz condition on <?j, we conclude by induction. 



< 



Remark 4.3 Assume that Xk is chosen to be an L 2 -optimal quantizer of Xt k for each k = 
0, ... ,m. It is in particular a stationary quantizer in the sense that E[At fc |A^] = Xk (see 
(15]). and by Jensen's inequality, we deduce that ||-Xfc|| 2 < ll^fc!^ - R- ecaum g (|4.15p . the 
inequality (|4.21j) in Theorem 14.21 gives 



max 



\Yl-Y k % < K(l + \\X \\ 2 )J2\\Xt 



X e 



l=k 



for all k = 0, . . . , m. In particular, if Xq = xq is deterministic, then Xq = xq, and we have 
an error estimation by quantization of the value function function for the discrete-time 
optimal switching problem at the initial date measured by: 

m 

max|«i(0,a;o)-flo(so)| < K{1 + \x \)^\\X tk - X k \\ 2 (4.24) 
l&q k=i 

Suppose that one has at hand a global stack of N points for the whole space-time grid, to 
be dispatched with N k points for each fcth-time step, i.e. X^fcLi ^k = N. Then, as in [2], in 
the case of uniformly elliptic diffusion with bounded Lipschitz coefficients b and a, one can 
optimize over the N k s by using the rate of convergence for the minimal L 2 -quantization 
error given by Zador's theorem: 



\x tk - x k \\ 2 



J24\Wk\\ 2 d 

j as iVfe — > oo, 



where (p k is the probability density function of Xt k , and \\(p\\ = (J \(p(u)\ r du) 7 . From [3], 

i 

we have the bound ||(/?fc|| 2 d < K^ft^, for some constant K depending only on b, a, T, d. 

d+2 

Substituting into (|4.24p with Zador's theorem, we obtain 



m 



tk 
l ■ 



max \vi(0,x ) - v l (x )\ < K(l + \x \, 
t&q k=l N k d 

For fixed h = T/m and N, the sum in the upper bound of the above inequality is minimized 
over the size of the grids T k , k = 1, . . . , m with 

d 

l 2(3+1) 



N, 



tt 



Em .2(3+1) 
k=l l k 
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where \x~\ := min{A;S N, k > x}, and we have a global rate of convergence given by: 

max \vi(0,x ) - v l {x ) I < _~ , 

h{Nh)d 

Actually even with no extra assumptions on b and <r, we have the same estimate, since 
for all r > 0, 

\\X tk -X k \\ 2 < C 2 , r \\X tk \\^ 1/d < KN~ 1/d , 
see Lemma 1 in |13| . 

By combining with the estimate (|3. 14|) . we obtain an error bound between the value func- 
tion of the continuous-time optimal switching problem and its approximation by marginal 
quantization of order hA when choosing a number of points by grid Nh of order l/h~. 
This has to be compared with the number of points N of lower order l/h d in the Marko- 
vian quantization approach, see Remark 14. 1[ The complexity of this marginal quantization 
algorithm is of order O (X^fc^=i ^k^k+i)- I n terms of h, if we take = Nh = l/h~ , we 
then need 0(l/h 3d+1 ) operations to compute the value function. Recall that the Marko- 
vian quantization method requires a complexity of higher order 0(l/h Ad+1 ), but provides 
in compensation an approximation of the value function in the whole space grid X. 



5 Numerical tests 



We test our quantization algorithms by comparison results with explicit formulae for op- 
timal switching problems derived from chapter 5 in [17j . The formulae are obtained for 
infinite horizon problems, that we adapt to our case by taking as the final gain the (dis- 
counted) value function for the infinite horizon problem. 

We consider a two-regime switching problem where the diffusion is independent of the 
regime and follows a geometric Brownian motion, i.e. b(x, i) = bx, a(x, i) = ax, and the 
switching costs are constant c(x,i,j) = Cij = 1,2. The profit functions are in the form 
fi(t, x) = e~P t kix" li , i = 1,2. From Theorem 5.3.5 in |17|). the value functions are given by: 



vi(Q,x) 
v 2 (0,x) 



A lX m + Kikix 
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X < Xy 



B 2 x m + K 2 k 2 x^ 2 - ci2, x>x\ 



A 2 x m+ +K 2 k 2 x^, 
A lX m+ + Kikxx'K - pal 
B 2 x m ~ +K 2 k 2 x^, 



X < X_2 

x_ 2 ^ X ^ x 2 , 

X > ~X* 2 



where A4, Bi, Ki, x^ and X2 depend explicitly on the parameters. In the sequel, we take 
for value of the parameters: 

b = 0, a = 1, coi = cio = 0.5, fci = 2, k 2 = 1, 71 = 1/3, 72 = 2/3, = 1. 

We compute the value function in regime 2 taken at Xo = 3.0 by means of the first 
algorithm (Markovian quantization). We take R = WXq and vary m, 5 and N. The results 
are compared with the exact value in Table 1. Notice that the algorithm seems to be quite 
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robust and provides good results even when 5m and ^ do not satisfy the constraints given 
by our theoretical estimates in Remark 14, 11 

In Table 2, we have computed the value with the marginal quantization algorithm. We 
make vary the number of time steps m and the total number of grid points N (dispatched 
between the different time steps as described in Remark I4.3f) . We have used optimal quan- 
tization of the Brownian motion, and the transition probabilities Trf were computed by 
Monte-Carlo simulations with 10 6 sample paths (for an analysis of the error induced by 
this Monte-Carlo approximation, see Section 4 in PQ). We have also indicated the time 
spent for these computations. Actually, almost all of this time comes from the Monte- 
Carlo computations, as the tree descent algorithm is very fast (less than Is for all the 
tested parameters). 

For the two methods, we look at the impact of the quantization number for each time 
step (resp. ./V and Nh) on the precision of the results. As our theoretical estimates showed 
(see Remarks 14.11 and 14. 3p . for the first method, increasing TV" higher than does not 
seem to improve the precision, whereas for the second method, we can see for several values 
of h that changing Nh from to h~ 2 or /i~ 3 improves the precision. 

Comparing the two tables, the first method seems to provide precise estimates with 
slightly faster computation times, and it has the further advantage of computing simul- 
taneously the value functions at any points of the space discretization grid X. However, 
since most of the time spent by our second algorithm was devoted to the calculation of 
the transition probabilities 7rju', if these were computed beforehand and stored offline, the 
marginal quantization method becomes more competitive. 



(m,l/8,N) 


62(0,3.0) 


Numerical error (%) 


Algorithm time (s) 


(10,10,10) 


2.1925 


3.0 


0.2 


(10,10,100) 


2.1863 


2.7 


0.5 


(10,10,1000) 


2.1852 


2.7 


1.4 


(10,100,1000) 


2.1882 


2.8 


8.5 


(10,100,5000) 


2.1882 


2.8 


40 


(100,10,100) 


2.1218 


0.31 


1.0 


(100,10,1000) 


2.1213 


0.33 


8.0 


(100,10,5000) 


2.1213 


0.33 


39 


(100,100,100) 


2.1250 


0.16 


8.6 


(100,100,1000) 


2.1250 


0.16 


82 


Exact value 


2.1285 







Table 1: Results obtained by Markovian quantization 
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(m,N) 


Y 2 


Numerical error (%) 


Algorithm time (s) 


(10,100) 


2.2080 


3.7 


4.4 


(10,1000) 


2.2174 


4.2 


4.9 


(10,10000) 


2.1276 


0.04 


5.8 


(100,1000) 


2.1233 


0.24 


36 


(100,10000) 


2.1316 


0.15 


48 


(100,50000) 


2.1301 


0.07 


65 


(1000,10000) 


2.1161 


0.58 


353 


(1000,50000) 


2.1213 


0.34 


498 



Table 2: Results obtained by marginal quantization 
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